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Superconducting States in Frustrating t-J Model: A Model Connecting High- Tc 
Cuprates, Organic Conductors and Naa;Co02 
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The two-dimensional t-J model on a frustrating lattice is studied using mean-field variational 
theories with Gutzwiller approximation. We find that a superconducting state with broken 
time-reversal symmetry (d-|-id state) is realized in the parameter region close to the triangular 
lattice. The frustration enlarges the region of superconductivity when t < for the hole-doped 
case, which is equivalent to t > for electron doping. We also discuss the SU(2) degeneracy at 
half-filling. The d-|-id state probably corresponds to the spin gap state at half-filling. 
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Frustration is an interesting problem and has been ex- 
tensively studied in spin systems. The effect of frus- 
tration in itinerant electron systems, however, is not so 
well understood. Resonating-valence bond (RVB) state, 
which is a liquid state consisting of spin-singlet bonds, 
was originally discussed in the Heisenberg model on a 
triangular lattice where the frustration was expected to 
destroy the long-range magnetic order leading to a quan- 
tum spin-liquid. It was speculated that the carrier 
doping to the spin system will induce the RVB states 
which are directly related to superconductivity.'^) Thus 
it is very important to investigate where the RVB states 
are realized in frustrated systems. 

Recently it has become possible to dope carriers in the 
two-dimensional organic conductor, k-(BEDT-TTF)2X, 
which has a lattice structure very similar to the trian- 
gular lattice. Actually it shows superconductivity under 
pressure.^' More recently, a superconductivity has been 
discovered in Naj;Co02 compounds, in which Co atoms 
form layered triangular lattices.^-' Although the super- 
conducting properties in these materials have not been 
clarified yet, it is very interesting to study the possibil- 
ity of RVB-typc superconductivity in a model which has 
the effect of frustration. In this paper we study the t-J 
model on a lattice as shown in Fig. 1, which interpolates 
between the square and triangular lattices. The organic 
conductors have such lattice structures. 

In the case of high-T^ superconductivity, hole doping 
to the square Heisenberg spin system induces (1^2 _y2- 
wave RVB-type superconductivity. ^"^"^ On the other 
hand, in the triangular lattice, the analysis of high- 
temperature expansion of the free energy suggests that 
the doping induces an RVB state when the hopping in- 
tegral t is positive for the hole-doped case.^^^ (Note that 
this corresponds to t < for electron-doped case.) When 
i < 0, the competition between the ferromagnetic ten- 
dency and singlet formation was observed. In this pa- 
per, in order to study the possibility of superconductiv- 
ity in the ground state, we use a variational theory with 
Gutzwiller approximation.^^ -'^^^ We find that the frus- 
tration enlarges the region of superconductivity when 



Fig. 1. Lattice structure of the frustrated two-dimensional t-J 
model. On the bonds forming the square lattice (solid lines), 
we have t and J terms, t' and J' terms are defined on the diag- 
onal bonds (dashed lines). 



i < 0, and that a superconducting state with broken 
time-reversal symmetry (chiral state) is realized in the 
parameter region close to the triangular lattice. We also 
discuss the SU(2) degeneracy at half-filling. 

The t-J model on a lattice in Fig. 1 is written as 

n = -tY, Pg{cIc,„ -f H.c.)Pg + JY.S,- s, 

-t'Y. Pciclcj, + H.c.)Pg + J'J2^^- '^J'(l) 

where Pq = Y\i{^ — ^i^nii) is the Gutzwiller projection 
operator which excludes double occupancies. and 
represent the summation over nearest-neighbor and 
one of the next-nearest neighbor pairs as shown in Fig. 
1, respectively. The sign of t is crucial since there is 
no electron-hole symmetry. Hereafter we use \t\ as 
an energy unit and study the parameter range, t' /t = 
J' / J = 0.0-1.0. The model becomes the square lattice t- 
J model when J' / J = 0; and the triangular lattce when 
J'/J=l. 

For this Hamiltonian, we consider Gutzwiller-type 
variational wave function Pg\^) with |<I>) representing a 
BCS mean-field wave function. For estimating the vari- 
ational energy, we use the Gutzwiller approximation in 
which the effect of the projection is replaced by a statis- 
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tical weight as follows^' 

(cjcj) = 9t{4cj)o, 
{Si ■ Sj) = gs{Si ■ Sj)o. 



(2) 



Here (• • • )o represents the expectation value in terms 

of 1$), and (• • • ) represents the norniaHzcd expectation 
value in In the simplest Gutzwiller approxima- 

tion we have 

_ 2(5 _ 4 

9t - iq:^' 9s - (i + j)2' (3) 

where 5 is the density of holes, i.e., the electron number 
is n = 1 — (5. It has been shown that this approxima- 
tion gives a fairly reliable estimation for the variational 
energy for the pure d2.2_j,2-wavc superconducting state 
when it is compared with the variational Monte Carlo 
results.^' 

In order to search for the state with the lowest vari- 
ational energy, we studied large unit cells assuming site 
dependent order parameters. In this Gutzwiller approx- 
imation, we can have order parameters. 



(4) 



We assume Xij^ = Xijl = Xij ^nd that Ajj is a singlet 
pairing, i.e., Aij = Aji. By minimizing the variational 
energy wc obtain a Bogoliubov-de Gennes equatiou"'^^^ 
and self-consistency equations. We have tried various 
trial states and found that the uniform state is always 
stabilized except for some special states realized only 
very near half-filling. Since these special states may be 
an artifact of the mean-field-type calculations, they 
are not discussed in this paper. The obtained phase di- 
agrams for J/t = 0.3, J' /t' = 0.3 are shown in Fig. 2 as 
a function of the ratio J'/ J and the electron density n. 

We find that in the region with J'/J < 0.675 the 
dj;2_y2-wave superconducting state is stabilized. This 
state is the same as in the square lattice case (J' = 0). 
When J'/J exceeds 0.675 a superconducting state with 
imaginary part of order parameter appears, which we call 
as d-|-id state as shown in Fig. 2. Since there is no contin- 
uous symmetry in our model, the notation of d-|-id-wave 
just represents a discretized version of da;2_y2-id,;,;-wave 
superconductivity. The detailed analysis is shown later. 
This state has excitation gaps all over the fermi surface 
and breaks time reversal symmetry. Another interesting 
feature is that for t < the superconducting region in- 
creases as J'/J becomes large, while its region decreases 
when t > {}. 

Before discussing the superconducting properties, let 
us here clarify the SU(2) degeneracy realized at half- 
filling in the present model. By solving the Bogoliubov- 
de Gennes equations we find various states at half-filling 
which have exactly the same variational energy with each 
other. For the uniform state, we have the energy disper- 
sion for quasiparticles as 




Fig. 2. Phase diagram obtained in Gutzwiller approximation of 
the RVB states in the frustrated two-dimensional t-J model with 
(a) positive t and (b) negative t. J/\t\ = 0.3 and J' /\t'\ = 0.3. 



where 

3 3 

£fe = - -^g s J iXx cos k^ + Xy cos ky) - -gsJ'XxyCos{kx + ky), 



(6) 

' 3 

' coskx + A^e'*^ cosky) + -gsJ'^xy coH{hj, + ky). 

(7) 



el + \Fk 



(5) 



Here A^, (A^) is the superconducting order parameter 
on the bonds in the a;-(?/-)direction, and A^.^, is in the 
diagonal direction. We denote {9') as the phase of the 
order parameter relative to that on the diagonal bonds. 
Note that, in this notation, the pure d2.2_j,2-wave state 
is represented as A^^ = Aj,, Aa;^ = and 9 = 9' = 1^/2. 

We find that, in each energetically degenerate state, 
these parameters satisfy the following relations; 

XxXv + Ax Ay cos(6i + 9') = 0, (8a) 

XxXxy + ^x^xy COS 9 = 0, (8b) 
XyXxy + \^xy COS 9' = 0. (8c) 

In this case the energy dispersion becomes 

Ek = \gs^Jj^{xl + A2) C0S2 fc, + J2(x2 + A2) COS^ ky + J'^{xly + ^ly 

(9) 
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The self-consistency equations are 
1 



1 



1 = 



1 ^ j9sJ cos^ 

gs J cos^ ky 
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N ^ 
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Ek 



1 jQsJ' cos^(fcx + ky) 
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Ek 



(10a) 
(10b) 
(10c) 
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Aj,, A^y, 6, 9') 



where N is the number of sites and the summation over k 
is in the Brillouin zone (in our model — tt < kx,ky < tt). 

Since there are 8 parameters {xx, Xyi Xxy, A 
and 6 equations (8) and (10), wc have remaining two de- 
grees of freedom which lead to the various states obtained 
at half-filling. This is an extension of the SU(2) degen- 
eracy discussed in the square lattice (J' = 0).8'i4) por 
example, one of the realized states at half-filling has 



D 

D' 



* 
* 
4 



0.5 



Xy = 0, Xxy 7^ 0, = Aj^ ^ 0, 



^xy 



TT 

4' 



0, 



(11) 



When Xxy — >■ 0, this state is naturally connected to the 
so-called d-|-is state realized in the square lattice. An- 
other state with 



Xx =Xy = = Ay 0, 

Xxy = 0, A^y ^ and 6* = 6*' = 



TT 

2' 



(12) 



is connected to the d3;2_j,2-wave state for the square lat- 
tice. 

Since these degenerate states have common values for 



D=^xl + Al = Jxl + A2, 



D' 



(13) 
(14) 



we plot their J'/ J dependence in Fig. 3(a). As J' / J 
approaches to 1, Z> becomes equal to D' . 

Next let us discuss the doped case. Once some holes 
are introduced, the above degeneracy is readily lifted, 
and the lowest non-degenerate state becomes the super- 
doncuting state which is shown in Fig. 2. We note that 
the flux phase proposed before'' is not stabilized. 
Since the order parameters are uniquely determined at 
finite 5, we determine the relative phase of the super- 
conducting order parameters, 9 and 6' , by taking the 
limit (5 — > 0. The obtained phase is shown in Fig. 3(b). 
When Axy starts to have a finite value at J'/ J = 0.675, 
the relative phases are 9 = 9' = 7r/2. This means that 
a small component A^y appears as an imaginary part 
(iAa;j^) in addition to the pure da;2_y2-wave state. As 
J' / J increases, the relative phase approaches to 27r/3, 
and when J' = J (triangular lattice) the three order pa- 
rameters Ax, Ay and Axy become symmetric with each 
other. This state is consistent with the other works. ^'^^ 
This gradual change of the relative phase gives a natural 
continuation between the d3;2_y2-wave superconductivity 
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Fig. 3. J'/J dependences of (a) order parameters and (b) relative 
phase = 0' , at half-filling. See details in the text. 



in the square lattice and 120-degree phase in the trian- 
gular lattice. 

In Fig. 4 we show the doping dependences of the or- 
der parameters A = A^, = Ay = Axy and x = Xx = 
Xy = Xxy determined self-consistently at /' = /. ,/' = J 
and J/t = 0.3 (triangular lattice). Qualitatively similar 
figure has been obtained recently, but the order pa- 
rameter X has not been taken into account, which plays 
an important role in the gauge theories for the t-J model. 
As shown in Fig. 4, the superconducting region is large, 
which is consistent with the high-temperature expansion 
study for the triangular t-J model when t> Q . Surpris- 
ingly the superconducting region is even larger for t < 0, 
in which case the high-temperature study has shown that 
the ferromagnetic instability takes place for J/t < 0.3.^°^ 
In order to clarify the actual ground state in the triangu- 
lar t-J model, more careful energy comparisons including 
ferromagnetism as well as the 120-degrees antiferromag- 
netism (realized at half-filling) are necessary. In partic- 
ular, there will be a competition between the ferromag- 
netic fiuctuation and singlet formations for J/t > 0.3 
which leads to a heavy fermion-like behavior appearing 
in the specific heat.^^^ The relation to this competition 
and the RVB states remains a future problem. In the 
region with very low electron density, we expect a triplet 
superconductivity. 

For the recently discovered superconductor, Na.TCo04, 
the sign of t is controversial.^®' If the sign is negative. 
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Fig. 4. Doping dependences of the order parameters in the trian- 
gular t-J model with (a) positive t and (b) negative t. J/\t\ = 
0.3. 



the superconductivity is more favored. For the organic 
superconductors, the sign of t is negative, but in this 
case, the electron density is just at half-fihing. In the 
t-J model calculation, we cannot discuss this supercon- 
ductivity and instead the Hubbard model has been used 
to explain the superconductivity via spin- fluctuation. 
In such calculations, only da;2_j,2-wave order parameter 
is observed. It will be an interesting problem to search 
for the d+id state in the parameter close to the triangu- 
lar lattice in the spin-fluctuation mechanism. 

There is another competition between antiferromag- 
netism and superconductivity in the phase diagram 
shown in Fig. 2. For the Heisenberg model, Weihong et 
al. discussed that the commensurate antiferromagnetism 
is stable for J' /J < 0.7 using series expansion meth- 
ods. This critical value is the same as found in Fig. 2, 
which means that the d-f-id state probably corresponds to 
the spin gap state at half-filling. It is necessary to com- 
pare the energies of antiferromagnetic state and d-wave 
superconductivity in the phase diagram. However it is 
well known that the mean-field theories overestimate the 
magnetic states. The Gutzwiller approximation docs not 
remedy this problem.^) In the square lattice, the varia- 
tional Monte Carlo studies have shown that the lowest 
variational state is a coexistent state of antiferromag- 



netism and d-wave superconductivity for 5 < 0.1.^^"^^^ 
As J'/ J increases and the magnetic frustration is in- 
duced, we expect that the region of antiferromagnetism 

shrinks and is confined in the very vicinity of half-filling. 
Therefore the RVB superconducting region in Fig. 2 will 
safely survive. Although the detailed study on this issue 
is necessary, the phase diagram shown in Figs. 2 and 4 
promises that an interesting time-reversal breaking su- 
perconducting state is one of the strong candidates for 
the ground state. 

In summary we have studied the frustrating t-J model 
within a Gutzwiller approximation. The frustration 
induces a superconductivity with broken time-reversal 
symmetry in the parameter region close to the triangular 
lattice. The relative phase between the order parameters 
gradually chang function of J' / J which connects 

the dj;2_j^2-wavc and typical 120-degree phase in the tri- 
angular lattice. Wc also discussed the SU(2) degeneracy 
at half-filling. Wc found that there is unexpectedly very 
large variety of states in strongly correlated frustrating 
systems. 

The author would like to thank very useful discussions 
and conversations with B. S. Shastry, G. Baskaran, and 
T. Koretsune. 
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